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Abstract 

It is well known that the linear mass spectrum of light mesons 
C**"* | in the large-iV c limit is dual to the perturbative QCD continuum. 

' We find the form of the linear spectrum which is maximally dual to 

^NJ ■ the perturbation theory. The obtained ansatz turns out to be the 

$— i ! spectrum of the Lovelace-Shapiro dual amplitude. This spectrum is 
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< 



chirally symmetric in the sense that it corresponds to zero values for 
the order parameters of chiral symmetry breaking in QCD. We fur- 
ther assume that the actual spectrum is dual to the maximally dual 
one within the QCD sum rules based on OPE and in this way recover 
(f) ■ some known important relations and statements. The experimental 

spectrum seems to support the assumed duality: The masses of reso- 
| nances cluster near certain equidistant values of energy squared and 

CsJ ■ the positions of clustering follow the maximally dual ansatz. The lat- 

ter fact signals the chiral symmetry restoration in the excited states. 
We discuss also two topics of recent interest: nonlinear corrections to 
' the string-like spectrum and the dimension-two gluon condensate. 

Ph! PACS: 12.90.+b, 12.38.-t,12.39.Mk 
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1 Introduction 



Thirty years ago A. Migdal [1] considered the problem of approximation of 
the perturbative asymptotics for the vector correlator (the so-called parton 
model logarithm) by an infinite sum of poles in the "best possible way". The 
obtained spectrum happened to be linear for masses, with the corrections to 
the perturbative logarithm being exponentially small. Quite unexpectedly, 
similar result was recently obtained in the framework of AdS/QCD [2]. At 
present many theoretical and phenomenological indications have been accu- 
mulated that with a good accuracy the spectrum is linear for masses square, 
i.e. it has a string-like form. In addition, according to the Operator Prod- 
uct Expansion (OPE) for two-point correlators [3], the corrections to the 
perturbative logarithm have to decrease polynomially rather than exponen- 
tially. On the other hand, saturating the two-point correlators by the linear 
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spectrum for masses square one can reproduce correctly the analytical struc- 
ture of OPE (polynomial corrections) with non-vanishing residues (constant 
in asymptotics). In this respect an interesting question arises: What linear 
spectrum for masses square approximates the asymptotics of correlators in 
the "best possible way"? Needless to say, the answer to this question is of 
importance for the effective string and AdS/QCD approaches. In the given 
paper we propose a solution of this problem within the QCD sum rules in the 
large- N c limit [4]. We show that in this case the properties of linear spectrum 
and OPE are completely determined by those of Bernoulli polynomials. This 
very fact allows to solve the problem. For the vector mesons the solution 
turns out to be nothing but the spectrum of the Lovelace-Shapiro (LS) dual 
model [5]. We demonstrate also that the deviations from this spectrum in 
the axial-vector channel can be ascribed to the spontaneous Chiral Symme- 
try Breaking (CSB). We propose a hypothesis that the actual spectrum is 
dual to the maximally dual one in the sense that OPE is identical for both 
cases. In a particular realization of this concept we rederive some important 
relations. The pattern of clustering of the experimental spectrum seems to 
be tightly related with the assumed duality. 

The paper is organized as follows. In Section 2 the general formalism 
of QCD sum rules in the large- N c limit is presented for the case of linear 
spectrum. In Section 3 we derive the maximally dual spectrum. Section 4 is 
devoted to the analysis of corrections to the linear spectrum. Section 5 con- 
tains some comments on the sum rules with finite number of states. Section 6 
deals with experimental spectroscopy. Dimension-two gluon condensate is 
discussed in Section 7. Some auxiliary material and discussions are given in 
Appendices. We conclude in Section 8. 

2 Sum rules for linear spectrum 

Let us briefly remind the formalism of QCD sum rules in the large- N c limit. 
In this limit the two-point correlation functions of quark currents are satu- 
rated by an infinite set of narrow meson states with the quantum numbers 
of these currents, i.e. they can be represented in Euclidean space as follows 



For instance, Ty = T7 M for the vector isovector case, where r are the Pauli 
matrices. On the other hand, their high-energy asymptotics are given by 




(1) 



2 



OPE [3] . We will consider only the vector (V) and axial- vector (A) isovector 
channels in the chiral limit. After appropriate subtractions one has 



..2 00 r 

fc=l 

The coefficients Ck only weakly depend on the scale /1 and momentum Q 2 
and can be calculated perturbatively. 
Define 

n^(g 2 ) = (-^g 2 + q,q v ) ti v > a (q 2 ). (3) 

The OPE reads as follows 

W > 12?r 2 Q 2 12ir Q* 9 Q 6 W7 ' 

where £ y = —7, £ A = 11. The symbols (G 2 ) and (gg) denote the gluon and 
quark condensate respectively. The residues are parametrized as follows 

Z VA (n) = 2F 2 A (n). (5) 

with Fv,a{ji) being electromagnetic decay constants (see [6] for details). 

The sum rules simply follow from comparison at each power of Q~ 2 of 
OPE (j3J) with the sum in Eq. ([T]) after summing up over resonances (in a 
chiral invariant way) and subtracting irrelevant infinite constants. 

Consider the linear ansatz for the mass spectrum (n — 0, 1, 2, . . . ) 



mj[n) 



an + m 2 , Fj(n) = F 2 = Const. (6) 



Such an ansatz is typically given by the strings, where the universal slope a 
is proportional to the string tension T: a = 2nT. Up to an infinite constant 
one has 

^ Q 2 + an + m 2 a V a J v 7 

n=0 v ' 

with ip being the Digamma function. Let us introduce the notations 

_ Q 2 _m 2 
z = — , x = — . (8) 
a a 

The ^-function has an asymptotic representation at large argument 

1 °° B 

lb(z + X] = In (Z + x) ; r— > — — tttt, (9) 

rv ' v ' 2(z + x) ^2k(z + x) 2k ' y ' 



Here B 2k denote the Bernoulli numbers (see Appendix A). Expanding Eq. © 
at z ^> x one arrives at (see Appendix B) 



*(„ + x) = In.- -ft^M, (10) 

fc=l 

where B k (x) are the Bernoulli polynomials. We could not find expansion (jlOp 
in the mathematical handbooks, so a derivation is given in Appendix B. 
This expansion is quite remarkable. It completely determines OPE for the 
linear spectrum through the Bernoulli polynomials and vice versa. Namely, 
equating expansions and ffTUl) one arrives at 

— l ^ + 2F E jfcQS + C0nSt= 12^V + ^' ( } 

One can always make logarithms equal by an appropriate choice of Const. 
Then 

^9 N c a 

= 24^' < 12 > 
and the general structure of sum rules for the linear spectrum is 

(-l) k N c a k B k (x) , . 

c *= (13 > 

Thus, the weak dependence on /x and Q 2 of the coefficients is replaced by 
a constant in the case of the linear spectrum. 

The sum rules considered in [6-15] are, in fact, particular cases of concise 
expression (jT3l) for k = 1, 2, 3 up to adding some separated resonances 0. We 
expect that one can trust these sum rules for k < 7. The reason is that the 
non-zero (i.e. even) Bernoulli numbers converge only up to B e , see Eq. (1601) . 
Then the divergence sets in reflecting the asymptotic nature of the expansion. 

Let us consider the vector case as an example. The first three sum 
rules (USD are (see Eq. (0}) 



= x-1/2, (14) 



a s (G 2 } 
YlTtF 2 a 
2e™ s (qq) 2 
3F 2 a 2 



x 2 -x + l/6, (15) 
x (x - 1/2) (x - 1) , (16) 



1 Similar expressions were derived also in [16] in a different way. 
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where F 2 is given by Eq. (fT2l . In the first sum rule, Eq. (|14[) . one has zero 
in the l.h.s. due to the absence of local gauge-invariant dim2 condensate in 
OPE. 

Consider the quantity in the l.h.s. of Eq. ( Tl6l) . with F 2 being substituted 
from Eq. (112)1 . 

lQ7r^ v a s (qq) 2 
e ~ N c a* ' 1 } 

(we remind that although a s and the quark condensate depend on the renor- 
malization scale, this dependence is almost negligible for the combination 
a s (qq) 2 ). Below we will be interested in the solution x = 1/2, corresponding 
to (qq) = 0. The quantity e shifts that solution, 

x=\ + 5. (18) 

Let us estimate how large is this shift due to a non-zero value of quark 
condensate. From the Gell-Mann-Oakes-Renner relation, 

m lfl = -(m u + m d )(qq), (19) 

with the standard inputs (the scale of the chiral symmetry breaking, 1 4- 1.2 
GeV is assumed) m u + m d « 12 MeV, f n « 92.4 MeV and m n « 135 MeV 
one obtains (qq) ~ —(235 MeV) 3 . The QCD coupling constant at the same 
scale is a s ~ 0.5. The last quantity we need is the slope a. If we accept 
an averaged value a « 1.14 GeV 2 from [37] then we obtain from Eq. ( fl7l) : 
e ~ 0.07. Substituting Eq. (1TBI) into the r.h.s. of Eq. ( fl6l) and neglecting the 
term C»(5 3 ), one arrives at 

5 ~ -4e. (20) 

It can be seen now that in the real world the corrections are of the order of 
5 w -0.3. 

Our subsequent analysis is valid, strictly speaking, for small corrections, 
\8\ 0.5, i.e. either in the limit of small enough quark condensate or in 
the limit of large enough slope. The latter means that the ground state is 
largely separated from the first radial excitation. Nevertheless, one can hope 
that the analysis is relevant for the real world: If one considers the planar 
limit of QCD (as we do), one regards the corrections 0(1/N C ) = (9(1/3) as 
small and this turns out to be a good approximation to the real world. In 
our case the magnitude of corrections is comparable with that of the planar 
limit. With this justification in mind, we will take i = 1/2 as a starting 
point in our analysis and will show that this works rather well for deriving 
various qualitative results. 
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3 Maximally dual spectrum 



It is well known that any linear spectrum is dual to the perturbation theory 
at one loop with free quarks since it reproduces the parton model logarithm 
in Eq. ffTTT) . However, in such a perturbation theory the power corrections 
are absent, while one cannot avoid them saturating correlators by the linear 
spectrum. In this situation it is interesting to find an ansatz which minimizes 
the power corrections. We will refer to such an ansatz as Maximally Dual 
Spectrum (MDS). 

First of all we note that from the absence of dim2 condensate, correspond- 
ing to C\ = in Eq. (jT5j) . follows the solution x = 1/2. Property (IMI) yields 
then C*2fe+i = 0. Thus, the half of power corrections is automatically zero. In 
particular, C3 ~ o s (qq) 2 = 0, i.e. one has no CSB. It is a good sign since the 
perturbation theory does not know about CSB. Relation ( 1631) shows, how- 
ever, that the zero values of C^+i condensates correspond to extremums of 
C*2fc ones and vice versa. Consequently, it is impossible to cancel all power 
corrections. In QCD this could mean that even if the chiral symmetry had 
not been spontaneously broken ((qq) = 0), some nonperturbative conden- 
sates would have existed. Presumably, these condensates are the purely 
gluon condensates with even powers of gluon fields, symbolically (G 2k ). We 
point out the negative sign of the gluon condensate provided by Eq. (TT3T) for 
x = 1/2, which indicates (if the quark-hadron duality is operative here) that 
the obtained spectrum cannot be the spectrum of real resonances, for mod- 
elling the latter one has to introduce some corrections. The value of gluon 
condensate is, however, very sensitive to such corrections. This fact can be 
used for constructing the corrections to the masses and residues, which are 
consistent with the physical values of OPE condensates [6]. In the present 
work we are interested in the ansatz, minimizing the power corrections in 
the OPE. We will not bother about the sign of condensates which it gives 
(as said above, it can be always improved by introducing small corrections). 
Of our concern will be a usefulness of this ansatz for a qualitative theoretical 
analysis. 

Thus, minimizing the half of power corrections one inevitably maximizes 
the other half. Let us formulate now the problem of minimization rigorously. 
For the time being it is convenient to measure the condensates in Eq. (|T3l) in 
the units of slope a, so ~ Bk ^ . Consider the first iV sum rules. We define 
the maximally dual spectrum as an ansatz providing the absolute minimum 
for the sum 




(21) 
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Consequently, one should solve the equation 



^- = 2^B k (x)B k . 1 (x), (22) 
fc=i 



and substitute the solutions into 

d 2 S{x) 



dx 2 



2 E ( B k-i(*) + ^B k (x)B k ^(x)^ 



(23) 



in order to to check the sign. Here relation (|63|) was exploited. The solution 
of the minimization problem turns out to depend on N. For N < 16 the 
minimum is delivered by x = 1/2. For N ^ 16 and uneven, the minimum 
is the same, but for the even N one obtains x ► It, f ), i.e. the roots 

N-+00 L 4 4 J 

of B 2k {x) in asymptotics (see Eq. (1671) ). However, as pointed out above, we 
cannot trust the sum rules at iV > 7 since the asymptotic divergence sets in. 
For this reason the last term gives the main contribution to the sums above 
since some N and, hence, completely defines their properties. Thereby we 
arrive at our corollary: The maximally dual vector linear spectrum is 

m 2 {n) = a (n + 1/2). (24) 

This is nothing but a LS-type of spectrum for the vector states. In fact, the 
same conclusion follows for any choice of minimality criterion of the kind 



S r (x) = 



V BJx) 



k=i 



k 



r > 0. (25) 



The maximally dual spectrum ( 1241) is chirally symmetric in the sense 
that the order parameters of CSB, (qq) and f n , are automatically zero due 
to Eqs. (fl4"|) and (fl6|) . This is not surprising since the perturbation theory is 
explicitly chirally symmetric. 



4 Corrections to spectrum 

The exactly linear spectrum is, of course, only approximation. In reality 
there are corrections as to masses, m 2 (n) = a(n + x) + A(n), as to residues 
F 2 (n). In the real world a smooth function A(n) hardly exists. Even in a 
much simpler case, dim2 QCD in the large- iV c limit, the 't Hooft model [17], 
the spectrum is obtained as a result of solution of the eigenvalue problem for 
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a highly nontrivial integral equation. As a result one has the linear spectrum 
only asymptotically, for the actual masses one observes seemingly random 
fluctuations around the asymptotic spectrum. In the real world one should 
have something similar. One can try to interpolate these deviations by a 
smooth function A(n). There are arguments [6] that for the self-consistency 
of sum rules these deviations should then decrease in n at least exponentially. 
An alternative way to demonstrate this point is presented in Appendix C. 
However, the self-consistency of sum rules admits also the existence of finite 
number of states, whose masses and residues are not subjected to a smooth 
in n parametrization. These states should be inserted "by hand" with the 
help of 5-function. This subject we consider in the present section. 
Each separated state contributes to the sum rules as follows 

*? _ F2 f^ (-l)W 

Q 2 + mf 1 Q 2 ^ 1 ) ' 

where index I enumerates these states. Each individual contribution is, gen- 
erally speaking, quite large in comparison with typical quantities appearing 
in the sum rules from the OPE side. On the other hand, the linear ansatz 
is a good approximation. Thus, in order not to spoil the sum rules it is 
reasonable to require that the sum over these contributions is approximately 
equal to the same sum over these states when they are subjected to the linear 
parametrization, i.e. we require the duality between the actual spectrum and 
the linear one in the sense of approximate coincidence of their expansions at 
large Euclidean momentum. This approximation we will denote by the sign 
" ~" . Let us separate in this way the first k states. One then has 

k k 

Y,F?mf^Y,F\l + x)\ 1 = 0,1,.... (27) 

1=0 1=0 

It is convenient to parametrize the deviations from the linear spectrum 
through the dilatations, 

if = XiF 2 , (28) 
m? = A t a(l + x). (29) 

When Xi = Ai = 1 one has no deviations. The requirement of duality ( 1271) 
reads then as follows 

k k 

MMT(i + x) { ^ J2( l + 1 = °> x ' • • • • ( 3 °) 

1=0 1=0 
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For the axial-vector mesons the first sum rule in Eq. (1301 has an additional 
pole at Q 2 = coming from PCAC due to propagation of the massless pion. 
This sum rule reads, 

k 

E A ^ fc > ( 31 ) 



p 



F 2 

1=0 

where f n = 87 MeV (in the chiral limit [18]) is the weak pion decay constant. 

Let us separate the first vector (p-meson) and axial-vector (ai-meson) 
mesons and assume that only these states can be strongly affected by CSB 
as their masses are close to Acsb- We have then for the sum rules with 
a = 0,1, 

1 ~ A p , (32) 

1 * J^ + K, (33) 
1 ~ A p>0l A Pj0l . (34) 

In addition, these relations are supplemented with Eq. (fi"2l) and the solution 
for intercept, x = 1/2, the latter means the duality between the real spectrum 
and MDS. 

In the real world with the chiral symmetry breaking there are some re- 
lations between the spectral parameters of the lowest states which indepen- 
dently hold in many realistic models and are in a good agreement with the 
experiment. In essence, they are consequence of "one resonance + contin- 
uum" ansatz. These relations are 

F 2 P « Vl (35) 
< ~ 2mJ, (36) 
24tt 2 
~N~, 



m l « — /?■ ( 37 ) 



The first relation, Eq. (l35il . is the KSFR relation [19]. The second one was 
originally obtained by S. Weinberg within the spectral sum rules [20] using 
Eq. (|35|) . Afterward it became clear that this formula expresses the fact of 
maximal mixing of longitudinal component of ai-meson with the pion after 
CSB, while between the vector and scalar states this does not occur [21]. 
Relation ( 1371) was originally obtained within the Borel sum rules [22]. At 
present it is understood as a consequence of vector meson dominance and 
usually it holds in models respecting this property (see, e.g., [23]). 

We also should respect relations (1351) - (1371) . The standard point of view on 
the sum rules is that they cannot completely determine the mass spectrum. 
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This is related to the fact that OPE is only asymptotic expansion for the 
correlators. The sum rules can give only some relations. To get predictions 
one should supplement them by some external inputs. First of all we observe 
that taking one of relations ( |35|) -(l371) as input, the other two relations are 
automatically reproduced by sum rules (|32l)-( l34|) . giving as a byproduct the 
relation for the slope, 

9 487T 2 „ 9 . , 

a ~ 2mJ ~ — (38) 

and for the axial- vector decay constant F ai ~ f n . For the parameters of 
dilatations this reads: \ p ~ A p ~ 1, A ai ~ 1/2 and A ai ~ 2. Thus, the mass 
and residue of p meson are not affected by CSB in our approximation. It is 
important to note that the obtained slope coincides with the one given by 
the LS model, Eq. i.e. we have not just a LS-type spectrum in Eq. (l2"4l) . 
but exactly this spectrum! 

However, in the axial-vector case the sum rules with % > 1 are now vio- 
lated. In particular, at i = 2 the contribution 2/^m* appears, giving rise to 
the condensate na s (qq) 2 ~ (4/j,-) 6 . In this way we reproduce the conclusion 
of Ref. [24]: For QCD in the large-N c limit spontaneous CSB with f^y^O nec- 
essarily implies the existence of non-zero local order parameters, where it was 
deduced from analysis of finite number of resonances within the sum rules. 
In our case this conclusion follows for an infinite number of states (which is 
more consistent with the large- N c limit), with important relations (!35|) - (l37|) 
being satisfied. 

In principle, one could add more separate resonances and obtain a better 
agreement with the experiment or accommodate the actual values of con- 
densates. However, the theoretical selfconsistency cannot be independently 
checked in this case because we do not know an analog of Eqs. fl35|) - fl37|) for, 
say, "2 resonances + continuum" ansatz. So we prefer to stop at this stage. 

We would like to point out that as a matter of fact we do not need any 
inputs (I35 i) -( 13~TI) . Due to PCAC they all can be reproduced automatically 
if we separate the ground axial-vector meson. Namely, if one neglects the 
exponentially small terms and possible contributions from higher orders of 
perturbation theory [25,26], the residues are given by the formula (see Ap- 
pendix C) 

where the factor c is determined from the coefficient in front of the pertur- 
bative logarithm, Eq. ffT2"]) . c = A" c /(247r 2 ). Let us deviate the lowest state 
from the linear spectrum, 

m 2 {n) = a [n + x + e5(n)] , (40) 
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where S(n) is the 5-function. Since the derivative of linear function is equal to 
the corresponding finite difference, one has from Eq. f[3"§]) for small variations 



F 2 (n) = c [m 2 (n + 1) - m 2 (n)] ~ ca [1 - eS(n)] . (41) 

Thus, even without sum rules one can see that if the mass of the ground state 
is increased in comparison with the value predicted by the linear spectrum, 
the decay constant of this particle is decreased, 

^ ~ - 1 - (42) 
5F*(0) ~ c [ ' 

Now let us write the same for the axial-vectors in terms of dilatations. The 
spectrum ( flOl) will be 

<(n) = a(n+^ + ^ll(5(n)). (43) 



Keeping in mind Eq. (|4T[) . it leads to 

A 01 ~ 1 + ^(1 -A ai ). (44) 

The system of Eqs. (il]) and (G3D has two solutions, {A 0l) A ai } = {1, 1}, {1/2, 2}. 
The latter is valid only if there is the additional contribution 1/2 to the first 
sum rule, Eq. f[3"3"j) . Identifying this contribution with the pion pole due to 
PCAC, all the results above follow immediately. Thus, the enhancement of 
the mass of axial- vector meson and the appearance of pion pole are intimately 
related. As a byproduct, Eq. f T44l) gives the constraint: m ai < \^3m p . 



5 Sum rules with finite number of states 

In this section we digress from the main line and consider the case of sum 
rules with finite number of resonances. There is a quite popular method of 
treating with such type of sum rules, the so-called Finite Energy Sum Rules 
(FESR) approach (for a review see, e.g., [27]). In short, the idea is to sepa- 
rate some (usually one) narrow resonances and make a cut-off at some energy, 
replacing the rest by the perturbative QCD continuum. These sum rules usu- 
ally work fairly well for differences of correlators, say U V (Q 2 ) — U A (Q 2 ), but 
for individual correlators they are not so good. A reason is that replacing 
the contribution of infinite number of states (any selfconsistent analysis of 
infinitely narrow mesons inevitably implies the large- N c limit, hence, the infi- 
nite number of states) by a cut-off is a quite rough interpolation. As a result, 
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one has an uncontrollable growth of condensates, C k ~ Yl F 2 (n)m 2 ^ k ~ l \n), 
which is unrealistic to compensate by the choice of cut-off at k > 2 in all 
sum rules simultaneously. In practice, the condensates, say C3 ~ a s (qq) 2 , 
are much smaller than even the first term, ~ F 2 (0)m 4 (0) in the given case. 
Consequently, the realistic condensates can appear only after some tricky 
cancellations of large numbers. Such predictions hardly can be stable in 
response to variations of parameters. 

We would like to propose another type of interpolation which is free from 
this drawback. Assuming linearity of spectrum (this is well settled exper- 
imentally, see, e.g., [28]) one can sum over resonances and arrive at sum 
rules f fl3|) . Then one can regard the r.h.s. of Eq. (Tl3]) as interpolating func- 
tions for contribution of finite number of states, fitting the parameters F 2 , a, 
and x. In particular, choosing a value for x near the roots of B 3 (x) naturally 
provides a small condensate C3. Such sum rules are automatically renor- 
malized during the summation procedure, as a result the dependence on the 
cut-off disappears. We expect that these sum rules should work at k > 2 in 
Eq. (fiHl) . when the complete set of roots of B k (x) on the interval [0, 1] sets 
in. 

Let us demonstrate an example of using such sum rules. Consider the 
third sum rule, 

C 3 (x) 1/2)0-1), (45) 

and assume the one- resonance saturation. Since in practice C^(x) is small in 
comparison with the typical values given by the r.h.s. of Eq. (f4"5l) (i.e. given 
by the hadronic part) this term can be introduced by small variations of x 
near the roots of the r.h.s. of Eq. ff^ol) : x — > x + e. Then 



<7 3 (l/2 + e)~~ C 3 ({1,0} + £)~|. (46) 



Thus 

W+j _i 

C 3 ({l,0} + e)~ 2' 
while from OPE, Eq. (T4j), one has 

ci = e = _^ 



(47) 



(48 



We do not bother about the exact absolute value in Eq. (j47j) since we have 
not varied F 2 and a. What is important for us is the sign in Eqs. (T4T|) 
and (1481) . The negative sign is possible only when x p ~ 1/2 and x ai ~ 1 (we 
do not consider the unphysical possibility x ai ~ 0). Therefore, we obtain the 
Weinberg relation, Eq. (1361) . making use of the third sum rule only. Thus, 
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the sum rules under consideration work in the region where the usual ones 
fail for the case of one separated resonance. Combining this relation with the 
Weinberg sum rules or FESR at k — 1,2 one arrives at the statement which 
is converse to that of cited in the previous section: For QCD in the large- 
N c limit spontaneous CSB with non-zero local order parameters necessarily 
implies the existence of f 2 7^ (a la Coleman- Witten theorem [29]). 

One can consider more resonances in Eq. (j45j) . i.e. regard the r.h.s. as 
an approximant for several states. This will result in the spectrum m 2 p (n) ~ 
a(n + 1/2) and m 2 ai (n) ~ a(n + 1), which is a particular spectrum of the gen- 
eralized LS amplitude, the so-called Ademollo-Veneziano- Weinberg (AVW) 
dual amplitude [30] . We remind that the LS dual amplitude is an amplitude 
for the reaction 7r + 7r — > 7r + 7r with linearly rising Regge trajectories and the 
correct chiral properties. It predicts the following spectrum for the p-meson 
trajectory and its daughters in the chiral limit 

m 2 LS (n, J) = 2m 2 p (n + J- 1/2), J > 0. (49) 

For J = one has to substitute J = 1 in Eq. ( ]49l . i.e. the scalar and vector 
mesons are degenerate. The AVW amplitude is a generalization of the LS one 
to the reactions 7r + A — > B + C. It contains also the spectrum for the pion 
trajectory (which includes the axial-vectors) and its daughters, predicting 

m 2 AVW (n,J)=2m 2 p (n + J). (50) 

The experiment does not support the AVW spectrum [31]: It fails in de- 
scribing the parity doubling on the daughter trajectories. Theoretically the 
reason is clear: The AVW spectrum leads to wrong chiral properties of corre- 
lators (the absence of chiral symmetry restoration at high energies, within the 
matching to effective models see [32]). Thus, this spectrum can be considered 
only as a low-energy approximation. For instance, it predicts a reasonable 
value for the chiral constant L 10 [18], 

L 10 = - 1 -[^(0)-U A (0)]=^ Xv) -' P{XA) 



2 



327T 

1.96 + 0.58 , . , 

4.4 -10" 3 , (51) 



32tt 2 

which is close to the prediction of one-resonance + continuum ansatz, 



1 



F 2 F 2 

P ai 



m 2 p m 2 ai j 



-2 + 0.5 , _ 



327T 2 



-4.7-10" 3 , (52) 



This happens because L\o is saturated at low energies. But if a physical 
quantity includes the high-energy asymptotics of correlators, for example, 
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the electromagnetic pion mass difference (for references see, e.g., [33]), 



Am, = - 3Q f2 dQ 2 Q 2 [U V (Q 2 ) - IT^Q 2 )] , (53) 

the AVW spectrum completely fails resulting in a divergence. 

Thus, sum rules ffT31 can be used for finite number of states, interpolating 
the hadronic contributions for k > 2. They supplement FESR at k = 1,2. 
Combining these two types of sum rules one can successfully describe the 
low-energy spectrum. 

6 Phenomenology and cluster duality 

The experimental spectrum of light nonstrange mesons is depicted on Fig. 1. 
The phenomenological analysis which led to that plot is performed in [34] 
(see also [35]). In short, the states below 1.9 GeV are taken from the Particle 
Data [36]. Above 1.9 GeV only a few states in [36] are cited. The resonances 
are taken from the Crystal Barrel Collaboration data [37]. The main reason 
is that this is the only experiment which performed a systematic detailed 
study of the energy region 1.9-2.4 GeV. For reliability, only the states seen 
at least in two reactions were selected. 

A prominent feature of the experimental spectrum is the existence of 
nearly equidistant clusters of states with a growing number of resonances. 
Qualitatively this coincides with the picture typically given by the strings, 
where the meson spectrum behaves as m 2 (n, J) ~ n + J (n is the number of 
radial excitation and J is the spin). Fixing the sum n + J = N one obtains 
the iV-th cluster containing N + 1 states. The chiral and axial symmetries of 
classical QCD Lagrangian, which are expected to be approximately restored 
at high energies (parity doubling), multiply the number of states in each 
cluster. 

It is convenient to introduce the cluster spectrum, M 2 (N) = A(N + X), 
which parametrizes the behavior of the experimental spectrum as a whole. 
The quantity A has the physical sense of mean slope and AX does of mean 
intercept. These parameters can be easily calculated from the experimental 
spectrum, they are presented in Table 1. We considered several cases. First 
of all, the vector and axial-vector isovector mesons, the sum rules for these 
channels are considered in the present work. Then all vector mesons, where 
the analysis is the same. Finally, all light nonstrange mesons. Moreover, 
as the first cluster contains only two states, it may be not well justified to 
consider this cluster on equal footing with the others, so we also displayed 
in Table 1 the data with the first cluster excluded. Since the data of Crystal 
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Clusters A, GeV 2 



X 



1+1- 
1+1- 
J = 1 
J = 1 
J = 1 



and 1 
and 1 
and 1 




0- 4 1.15 

1- 4 1.14 
0-2 1.11 

0- 4 1.12 

1- 4 1.13 
0-2 1.10 

0- 4 1.13 

1- 4 1.13 
0-2 1.14 



0.54 
0.55 
0.59 
0.53 
0.51 
0.56 
0.54 
0.55 
0.54 



All mesons 
All mesons 
All mesons 



Table 1: Parameters of cluster spectrum m 2 (N) = A(N+ A), N = 0, 1, 2, 3, 4 
in Fig. (CQ). 

Barrel (the last two clusters) are still absent in Particle Data, we show also 
the results without this data (the first three clusters). It turns out that the 
results vary only a little. 

A striking feature of experimental spectrum is that on average it behaves 
as MDS, with the slope being very close to the predicted one, Eq. (13*51) . Since 
MDS is chirally symmetric, we conclude that although CSB shifts significantly 
the masses of some states, when one considers the whole light nonstrange me- 
son spectrum with more than 80 resonances, i.e. when the statistical weight of 
these several states is small, one observes the approximately chirally symmet- 
ric spectrum, as if CSB were absent at all! The experimental spectrum also 
reveals the universality of linear spectrum for any quantum numbers, which 
seems to allow the generalization of the conclusions for vector channels made 
in this paper to arbitrary ones. It must be emphasized that Particle Data 
(the first three clusters) contains enough states to arrive at the conclusion 
above. The Crystal Barrel Data convincingly confirms this observation. 

How the existence of approximately equidistant clusters is related to the 
hypothesis of duality of real spectrum to MDS? If deviations from MDS are 
sporadic to both sides, then the duality naturally leads to the clusters with 
the positions determined by MDS. However, we should take into account the 
possibility of systematic deviations from MDS to one particular side, which 
should be compensated by deviations to opposite side in other channels in 
order to provide a stable cluster. This can be achieved by generalization of 
duality considered above in the following way 




(54) 
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where I enumerates all possible sets of quantum numbers with fixed spin J. 
The spin defines a Lorentz tensor structure in front of the correlator, like 
in Eq. ([3]) for J = 1. This tensor structure is extracted in Eq. ( jMJ) - The 
mean correlator in the l.h.s. of Eq. (151]) is supposed to be saturated by 
MDS, hence, defining the positions of clustering the poles in the r.h.s. This 
generalized duality we would call the cluster duality. The case considered 
in Section IV is a particular case of cluster duality for two channels, when 
the spectrum of one of them (vector mesons) is locally dual to MDS (i.e. 
coincides). Then Eq. ( 1541) requires at least the global duality to MDS for the 
axial-vector states. 



7 Dimension-two condensate 

The authors of Ref. [38] advocated the idea that although the local dim2 
gauge invariant operator is absent in the standard OPE in the chiral limit, 
one can construct a gauge non-invariant dim2 operator from the gluon fields 
and the minimal value of its vacuum expectation value (achieved in the Lan- 
dau gauge), A 2 ~ (G), could have a definite physical sense, encoding some 
important nonperturbative effects. This subject soon got a certain popu- 
larity (see, e.g., references in [15]). Usually one attempts to improve the 
agreement with phenomenology introducing this condensate and in such a 
way one obtains an estimate on its value. In particular, within the sum rules 
in the large- N c limit some estimates have been derived in [12, 15] (in [6] it 
was noticed, however, that within this theoretical setting the dim2 conden- 
sate hardly can be detected). Let us pursue this way, i.e. let us assume the 
existence of dim2 gluon condensate and try to estimate its value within the 
linear spectrum. 

Consider the vector channel. The dim2 gluon condensate modifies the 
first sum rule [39], Eq. dTJ) 

a A 2 

-^^2F\x-l/2). (55) 

The axial- vector channel will have the same contribution, but the appearance 
of / 2 hampers the estimation in that case. As we have shown (see Table 1), 
in reality Ax = x — 1/2 > 0. Since phenomenologically A 2 < ("effective 
tachyonic gluon mass") we have agreement in the sign. Suppose that the 
deviation from the maximally dual ansatz, x — 1/2, happens completely due 
to existence of A 2 7^ (this yields disagreement in the sign for the quark 
condensate term in the third sum rule, Eq. ffTB"]) . but for the time being we 
do not bother about that). We expect that this should give a significant 
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overestimation for the quantity in question. Using Eq. (|12p one arrives at 

a s X 2 N c 

~ — aAx. (56) 

7T 3 

Substituting the approximate experimental values for x and a from Table 1 
into Eq. (I56p one obtains the estimate — ~ 0.03 -h 0.06 GeV 2 (we neglect 
the extremal values). This supposedly overestimated value turns out to be 
less up to one order of magnitude than ones, usually cited in the literature. 
In fact, the experimental deviations Ax have the same smallness as the chiral 
corrections. Say, in the LS amplitude the pion mass yields the correction 

x ts = -+ 2 77 2 « 0.53. (57) 
2 m l p — m£ 

Thus, in the chiral limit we cannot make the case in estimating quantitatively 
the dim2 condensate within our theoretical setting, let alone the fact that 
the chiral corrections are usually less than 0(1/N C ) ones. As a matter of 
fact, our estimate is well consistent with the absence of dim2 condensate. 



8 Conclusions 

We give a new insight into the phenomenological success of the Lovelace- 
Shapiro dual amplitude: Its spectrum is maximally dual to the perturbation 
theory among the linear spectra. This spectrum turns out to describe the 
positions of clusters of meson states near some values of energy. This indi- 
cates on the existence of certain duality between the real spectrum and the 
maximally dual one. The reason seems to be the chiral symmetry restoration 
in the excited states since the maximally dual spectrum corresponds to zero 
order parameters of the chiral symmetry breaking in QCD. In the given pa- 
per we have demonstrated how this duality may work within the sum rules 
giving quite encouraging results. Our analysis shows that there is a hope to 
describe quantitatively the global properties of observed light meson spec- 
trum within the QCD sum rules, at least with the accuracy typical for the 
large- N c limit. 

We would like to mention also that the QCD sum rules suggest a possible 
way to prove the asymptotic linearity of the light meson spectrum, at least 
for the vector channels. Namely, due to relation ( 1771) this is tantamount to 
an independent demonstration that the decay constants asymptotically tend 
to a constant value. 

A possible continuation of the present work is an analysis of channels with 
arbitrary spins. The experimental spectrum shows that the results must be 
similar. 
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Appendix A. Bernoulli numbers and polyno- 
mials 

The Bernoulli numbers can be represented as follows 

= (-^"fpC^), (58) 
B 2k +i = 0, k>0. (59) 
Here ( is the Riemann function. The first few numbers are 

B = 1, B 1 — --, B 2 = -, 

The Bernoulli polynomials are defined as 

B k (x) = J2(! n ) B ^ m - (61) 

m=0 ^ ' 

The first few polynomials are 

B (x) = l, B 1 (x)=x-\, B 2 {x) = x 2 - x + \, 
3 1 

B 3 (x) = x 3 - -x 2 + -x = x(x - 1/2) (x - 1), 

Zi z 

B 4 (x) = x 4 - 2x 3 + x 2 - — . (62) 

oU 

The Bernoulli polynomials have many interesting properties. Below we cite 
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those which are relevant for our analysis. For k > one has 



B' k {x) = 






(63) 


B 2k+1 (0, 1/2, 1) = 


o, 




(64) 


B 2k (0, 1) = 


E>2k, 




(65) 


a»(i/2) - 


> — 

->oo 




(66) 


fe- 


>~ 

->oo 


cos(27nr), 


(67) 


5 2fc +i(x) 

fe- 


>~ 

->oo 


sin(27ra;). 


(68) 



It should be noted that in fact asymptotic properties (I6"6"l) - fl6"8l set in quite 
rapidly Say B 2 {\) = -|-B 2 , but already B±(\) » -0.975 4 - As to Eq. (1671), 
we are interested in the roots on the interval [0, 1] only. Already for B2(x) 
the approximate roots 0.21 and 0.79 are close to the asymptotic ones, \ and 

3 
4' 



Appendix B. Expansion of i/> function in Bernoulli 
polynomials 

In this Appendix we derive Eq. (ITU]) . Using property (15*9"]) one can cast Eq. © 
into a form 

-ib(z + x) = -lnz-ln (l + -) + -— ^ — - + V - . Bk . . , (69) 
^ v ; V zJ 2{z + x) j^k{z + x) k ' K 1 



Expand the last three terms in the r.h.s. of Eq. (1691) at z ^> x, 

fe=i fc=i v 7 



2(z + x) L ) z k 2^ kz k \i niX 

v ; k=0 k=l v 7 

k(z + x) k ~ ^kz k ^ '\ I J z l 

k=2 K ' k=2 1=0 v 7 

^ 1 * (-l) k ~ m (k- 1 



z k ^— ' m \k — m 

k=2 m=2 



k—m 



k=2 m=2 
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(71) 



In the last case we multiplied the sums according to the formula 

oo oo oo n 

E a ^ n E b k( n )y k = E yU E °A-fcW, 

n=2 k=0 n=2 k=2 

then used identity ( k Z^) — ir( m ) anc ^ relation ( J59l) . Summing Eqs. ( I7U1) . 
(EE]) , ( 1721) and substituting the result into Eq. (1691) one obtains finally 

OO / fc / / \ 

- ^ + = - In z + £ ^ E J ^^" m - ( 73 ) 

k=l m=0 ^ ' 

Using definition of Bernoulli polynomials ( 16TI) one arrives at Eq. ( jTUi) . 



Appendix C. Systematic nonlinear corrections 
to spectrum 

Let us analyze the constraints on the possible nonlinear corrections to the 
string-like spectrum. A naive way is to expand in powers of corrections in 
the sum of Eq. (J7J). However, the validity of this procedure requires a strong 
convergence properties and apriori we do not know them. Instead we will 
make use of the Euler-Maklaurin summation formula, which is an asymptotic 
series for a sum 

J2f(n) = / f(x)dx-B 1 [f(0)+f(N)} 
Jo 



n=0 



+ Eff [/^(^-/^(O)]. (74) 

k=2 

We are interested in functions with the property f( k \N) > 0, k 

0, 1, . . . , for which Eq. ( 1741) reads 

-^- 1} (0). (75) 



N^oo 



J " poo _ _ R 

E/w = / /(^-Ett/ ( 

n=0 ^ k=l 



In our case 

/(») TTT . (76) 

v ' z + n + x + A(n) v 7 

For the linear spectrum, F 2 (n) = const, A(n) = 0, expression (1751) gives 
the asymptotic expansion of Digamma function, Eq. (Q. In particular, the 



20 



logarithm appears from the integral. This imposes a constraint on F 2 (n): 
To reproduce the analytical structure of OPE (i.e. to have the series in z~ k 
only) the residues has to behave as 

F 2 (n) ~ ^M, (77) 

up to some exponentially small terms [6]. We will not be interested in the 
latter. The problem is reduced to the corrections to the linear mass spectrum, 

/(B ) „ 1 + A, '») - (78) 
v ' z + n + x + A(n) v ' 

Substituting Eq. ( 1751) into Eq. ( 1751) and expanding in inverse powers of z one 
obtains an infinite series at each power because now A^(0) 7^ at k > 1. 
Say, at 1/z one has a contribution 

1 00 R 

i: ~£*A<M>(0). (79) 

fc=l 

Now we require that small corrections do not have to generate uncontrollable 
divergences in all results. If we want to preserve the present scheme of 
analytical treatment with the sum rules we need to impose the convergence 
of the sums like (1791) . From relation (1581) for large even k one has then the 
following condition in the case of sum ( 1791) 

A«(0) 



1 

2^ 



A( fc - 1 )(0) v ; 

For power-like corrections, A(n) ~ n~ a (a > 0), condition (IHUj) cannot be 
fulfilled, as well as for the corrections like A(n) ~ n a e~ l3n {(3 > 0), where 
a is not a positive integer. But for purely exponential ones, 5(n) ~ e _/3n , 
condition (IHUj) can be satisfied if /? < 2tt. The same is true for a special type 
of "mixed" corrections: A(n) ~ n a e~^ n (a, (3 > 0), where a is integer. At 
other powers of 1/z the situation is more involved. Let us give only estimate 
for the exponential corrections (they have a property A( fc )(0) = (A'(0)) fe ), 

^M^! aM(0)i (81 ) 

fc=l 

where c m = 0(1), c m > 1. One obtains then (3 < 2n/c m for convergence. 

Thus, within the present scheme there are strong constraints on the 
derivatives of corrections. The power-like corrections cannot satisfy them, 
but the exponential ones do can. On the other hand, the exponent (3 is 
also constrained. Actually, it can be so small that in the beginning of spec- 
trum the decreasing is slower than, say, 1/n. The presented estimates are in 
agreement with the results and fits obtained in [6]. 
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Figure 1: The spectrum of light non-strange mesons from [36] and [37] (for the 
last two clusters) in units of m^ 770 y Experimental errors are indicated. Circles 
stay when errors are negligible. The dashed lines mark the mean (mass) 2 in each 
cluster. The clustering occurs at 0.78, 1.32, 1.70, 2.00, and 2.27 GeV. The following 
states are displayed (in MeV): vr: 140, 1300±100, 1812±14, 2070±35, 2360±25; 
fa 1200-1500, 1770 ±12, 2020 ±38, 2337 ±14; r?: 1294 ±4, 1760 ±11, 2010^, 
2285±20; a : 1474±19, 2025±30; p: 775.8±0.5, 1465±25, 1720±20, 2000±30, 
2110 ± 35, 2265 ± 40; a x : 1230 ± 40, 1647 ± 22; uj: 782.59 ± 0.11, 1400 - 1450, 
1670 ± 30, 1960 ± 25; fv 1281.8 ± 0.6, 1971 ± 15, 2310 ± 60; h x : 1170 ± 20, 
1965 ± 45, 2215 ± 40; h: 1229.5 ± 3.2, 1960 ± 35, 2240 ± 35; m: 1376 ± 17, 
2013 ± 25; f 2 : 1275 ± 1, 1638 ± 6, 1934 ± 20, 2001 ± 10, 2240 ± 15, 
vr 2 : 1672 ±3, 2005 ±15, 2245 ±60; r/ 2 : 1617±5, 2030 ±16, 2267±14; 
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1653_ 15 , 
2293±13; 

a 2 : 1318.3 ± 0.6, 1732 ± 16, 2030 ± 20, 2255 ± 20; 
1975 ±20, 2195 ±30; f 3 : 2048 ± 8, 2303 ± 15; uj 3 : 
2285 ± 60; p 3 : 1688 ± 2.1, 1982 ± 14, 2260 ± 20; 
2025 ± 20, 2275 ± 25; 63: 2032 ± 12, 2245 ± 50; 
2283±17; p 4 : 2230±25; a 4 : 2005^, 2255 ±40; 
p 5 : 2300 ±45. 25 



p 2 : 1940 ±40, 2225 ±35; cj 2 : 
1667±4, 1945±20, 2255±15, 
a 3 : 2031 ± 12, 2275 ± 35; h 3 : 
7r 4 : 2250 ± 15; / 4 : 2018 ± 6, 
774: 2328±38; w 5 : 2250±70; 



